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PACS. 36.40. -c - Atomic and molecular clusters. 
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Abstract. - The periods of the harmonic oscillations of the ion core of charged sodium clusters 
around the equilibrium shapes are considered. It is found that these periods are of the order 
of magnitude of the experimentally measured relaxation times of the plasmons, which suggests 
the importance of the electron-ion coupling and stresses the role played by the electron-phonon 
interaction in the dissipation of the plasmon energy. The relation of the process to fission is 
briefly discussed. 



Introduction. - A considerable progress in cluster physics has been achieved for the past 
years. The picture of the collective motion of the electron system is well understood, [1-3]. 
The collective electron modes have also been studied experimentally. On the other hand, 
experiment shows a strong electron-ion correlation, with a characteristic time of ^ 10^^^ 
s [4]. There is no complete picture to explain such a strong correlation occurring in spite of 
the huge difference in the masses which actually justifies the use of the Born-Oppcnheimer 
approximation. In principle, the energy of a plasmon can be directly transferred to the phonons 
within a time which follows from the analytic solution of a simple model to be presented in 
greater detail elsewhere, 

where ujpi is the plasmon frequency, M is the mass of the ion and m is the electron mass. 
For sodium clusters, this time is of the order of 10~"'^^s, which may be even an order of 
magnitude shorter than the experiments show. The relaxation time of the plasmon may be 
identified with Tpi, eq.Q). We summarize now the derivation of eq. Q). We consider plane 
waves propagating along the z axis. We denote the displacements of the electronic and ionic 
distributions, respectively, hy u = u{z,t) and v = v{z,t). We denote the fluctuations of 
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the electronic and ionic charge densities, respectively, by 5pe = and 5pi = — e|^. The 
Lagrangian describing the plasmon-lattice dynamics reads 




where u' = u{z'), v' = v{z'), m, is the electron mass, M is the ion mass and ^oe{^)'^ is the 
elastic energy density of the lattice. The equations of motion read 

-M— + a— + Aire^u -v) = Q 

ot^ oz^ 

-TO— +47re2(z;- m) = 0. (2) 

It follows that, for an appropriate value of the wave vector, the number of plasmons at time 
t is 

TVp; = uu* cx cos^ (^y^^P'*) ' 

so that the plasmon energy will be completely transferred to the lattice in time Tpj. 

Our present purpose is to draw attention to the fact that the periods of the collective 
vibrations of the ionic core remarkably coincide with the detected electron relaxation times. 
This fact supports the assumption that the interaction of the electron system with the collec- 
tive ionic modes plays a rather essential part in dissipation, as, e.g., in the electron-phonon 
interaction in crystals. There is no need to assume that the phonon excitation will necessarily 
result into fission, though the two aspects are generically related to one another, so that fission 
is the ultimate form of a superposition of many of phonons in the limit of large amplitude 
oscillation [5]. The importance of the collective modes for the electron relaxation was also 
noted in ref. [6]. The influence of the ionic degrees of freedom on the electronic excitations 
has been considered in [2,3,7]. 

In turn, cluster fission is a process of great interest. Experiment shows predominance of the 
strongly asymmetric fission accompanied by emission of mono-, di- or, more rarely, trimers. 
Symmetric fission remains among the most important topics of research. Study of symmetric 
fission allows one to better understand the dynamics of the interplay between single-particle 
and collective degrees of freedom [5] . 

In paper [8], we have considered the interplay between rotational and vibrational modes 
of the collective motion in clusters, finally leading to fission. It was specifically found that 
rotation favors fission through phase transitions occurring in the shape of the clusters rotating 
with large angular momentum. This is similar to nuclear fission. On the other hand, only 
neutral clusters were studied in [8] . It is well known, however, that the charge of the clusters 
plays an important role in fission. 

In the present paper we extend our considerations to charged clusters. Moreover, we 
compare the obtained results with the experimental data, which became available after paper 
[8] was submitted for publication. Sodium clusters with N = 18, 43, 92 and 470 atoms are 
considered. In the light of recent experiments [4], we discuss the role played by the collective 
motion of the core in the electron-ion c;orrc;lations. Onr discussion is based on the LDM (sec [9] 
and references cited in [5]) which is enough for the present purposes. We take into account 
surface and Coulomb energies. The inclusion of a shell-correction term generally allows one 
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to obtain a detailed description of specific features of clusters, such as binding energy per 
particle, ionization energy etc. [5]. Here, we leave out these aspects, focusing rather on 
general tendencies than on detailed descriptions of particular properties of individual clusters. 
On the other hand, the LDM is well-suited for the description of fission, due to its intuitive 
appeal and transparency. Many properties of nuclear fission have been understood in the 
framework of the LDM. Representative sodium clusters in a large range of the number N of 
the constitutive atoms is dealt with, one of which, N — 92, is close to the clusters studied 
experimentally in [4]. We also draw special attention to possible cases of soft clusters with 
charge close to the critical value, just at the border of Coulomb stability, namely, clusters 
with N = 43, q = 3 and N = 470, q = 10, q being the cluster charge. 

Outline of the model. ~ We consider oscillations of an incompressible liquid drop such 
that, at each instant, the shape is that of an axially-symmetric spheroid, with time-dependent 
half- axis c in the direction z, and time-dependent half- axes a — b in the perpendicular direc- 
tions X and y. We assume irrotational flow. In view of volume conservation, the values of c 
and b are related, at each instant, by cb^ — , R being the radius of the equilibrium spherical 
shape. For sodium, R = 3.93iV^/'^bohr, N being the number of the atoms in the cluster [12]. 
We choose c as the collective variable. The mass parameter can be found by solving the 
Laplace equation for the velocity field with appropriate boundary condition, analogously to 
the case of small-amplitude multi-pole vibration [10]: 



(4) 



The velocity field is 

-v{x,y,z) = -Vxix,y,z) . 
A solution of eq. satisfying the proper boundary condition is 



(5) 



X 



4c 



(6) 



(7) 



The vibrational kinetic energy can be found as follows: 

T = ^m/ (Vx)'dl/, 

where denotes the mass per unit volume. From Eq. © one gets the expression 

T=^Mic)c\, M{c)^^Mil + ^u') , (8) 

where M is the total mass of the cluster, and u stands for [R/c). 

The potential energy of deformation arises from the interplay of the opposite effects of the 
restoring surface and the repulsive Coulomb energies. For the surface energy, the following 
expression was derived in [8] , 



Vsurf{.c) = 2wi?' 



V{c)surf = 27raR^ 



^1/2 arcsin vl-^? 



arcshV u-^ — 1 



- inaR^ for R < c , 



AnaR"^ for R> c , 



(9) 
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where a is the surface tension. For sodium clusters, we take a = 3.8 x 10^'^ eV / bohr^ in 
agreement with the Stabihzed JeUium Model [11,12]. The capacities of a prolate spheroid with 
half axes c> a — b and of an oblate spheroid with half axes a = b > c are [13], respectively, 



and C^^^^. (10) 
cosh-^(c/a) cos-i(c/a) ^ ^ 

These formulae allow us to obtain the Coulomb energy of a charged spheroid, 

Vcoulic) = ^. (11) 

For small deformations the potential energy becomes 

V^'Hc) -K% + V^t, (12) 

Vsulfic) =\hur}{c-Rf (13) 

1. 



^CoLW ^-i^kcoui{c~R)\ (14) 
where terms of higher order than the second in (c — R) have been neglected and 

ksurf = = 0.0382el//a2 , kcoui = = 0-1794ey-|^ . (15) 

flop's ^^'^O 

Here, is the Bohr radius, = 3.93, the mean radius per particle being rgUo, and q = 
0, 1, 2, • • • is the charge of the cluster. As in [8], the kinetic energy of the oscillating cluster is 

^-ix(0)c^ (16) 

where A^(0) = 6.417 • lO^eVN and the oscillation period is 



^ = 2^1/ 1 ^ • (1^) 

For small amplitudes, our assumptions result in the quadrupole vibrations of a liquid drop. If 
these assumptions are used in conjunction with different velocity fields from the one considered 
in the present paper, other types of fission will be favored, as, for instance, asymmetric fission. 

Results and discussion. - The results of the calculation are presented in Table for 
typical cluster with numbers of atoms N — 18, 43, 92 and 470, the charge q being such that 
the clusters stay within the stability limit. 

The periods can be compared with the experimental data [4]. In that paper, characteristic 
times for the relaxation of the plasmon energy were obtained, which turned out to be 2.5 x 
10^^^ s. As one can see from Tabled the plasmon relaxation time is about the same as the 
characteristic period of the collective vibration. Naturally, for the same period, the plasmon 
makes thousands of oscillations. The plasmon collective energy is quite likely dissipated by 
inducing the excitation of the cluster phonons. The electrons collective energy is directly 
transferred to the ion core, exciting phonon degrees of freedom. In view of the relation 
ijjph ^ ^ph the number of excited phonons is very large, n ^ 1 [5]. Such many-phonon 
coherent excitations can then serve as doorway states for fission [15], unless they are destroyed 
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Table I - Behavior of sodium clusters with 18, 43, 92 and 470 atoms. In the second column, the 
charge q of the cluster is shown. For each cluster type, the constants ksurf o,nd kcoui, characterizing, 
respectively, the surface and the Coulomb energies, and the period t of small amplitude oscillations 
around stable configurations, are presented. 



Species 


q 


ksurf (eVag ) 


kcoul (eVa^ ) 


7, 7 

^surf rZCoul 


r ( s) 


Nai8 





0.382E-01 





0.382E-01 


0.193E-11 




1 


0.382E-01 


0.997E-02 


0.282E-01 


0.224E-11 







0.382E-01 





0.382E-01 


0.298E-11 


Na43 


1 


0.382E-01 


0.417E-02 


0.340E-01 


0.316E-11 




2 


0.382E-01 


0.167E-01 


0.215E-01 


0.397E-11 




Q 
O 


U.ooZrj-Ul 


U.o (O-tj-Ui 


U.OOirj-Uo 









0.382E-01 





0.382E-01 


0.436E-11 




1 


0.382E-01 


0.195E-02 


0.362E-01 


0.448E-11 


Na92 


2 


0.382E-01 


0.780E-02 


0.304E-01 


0.489E-11 




3 


0.382E-01 


0.175E-01 


0.206E-01 


0.593E-11 




4 


0.382E-01 


0.312E-01 


0.700E-02 


0.102E-10 







0.382E-01 





0.382E-01 


0.985E-11 




3 


0.382E-01 


0.344E-02 


0.348E-01 


0.103E-10 




5 


0.382E-01 


0.954E-02 


0.287E-01 


0.114E-10 


Na470 


7 


0.382E-01 


0.187E-01 


0.195E-01 


0.138E-10 




8 


0.382E-01 


0.244E-01 


0.138E-01 


0.164E-10 




9 


0.382E-01 


0.309E-01 


0.728E-02 


0.226E-10 




10 


0.382E-01 


0.382E-01 


0.298E-04 


0.353E-09 



by dissipation. The latter process dominates if the spreading width Tg is greater than the 
mean spacing of the n-phonon states, fu^ph, 

Ts^hujph . (18) 

In the case of nuclear fission, the opposite relation to H18(l . namely, tiujpin takes place [15]. 
In this case, the n-phonon structure is experimentally observable. This corresponds to the 
Frenkel picture of fission [14], where the collective deformation arises from a large amplitude 
oscillation, which, due to an-harmonic effects, never comes back to the initial position, but 
proceeds towards fission. In this picture, fission occurs on a competitive basis with dissipation 
of the phonon collective energy inside the ion core. The observed fact that fission most 
likely occurs through the very asymmetric mode of evaporation of monomers or dimers means 
that strong dissipation of the core collective motion takes place, and that the relation H18|l 
is satisfied. Moreover, from the expected phonon dissipation time Tdiss ^ 10~^^s, we can 
estimate the related spreading width of the n-phonon states which turns out to be 

r^' ~ 7 X IQ-'^eV . (19) 

That is, the many-phonon states strongly overlap. This means that classical approaches, such 
as the present one, are appropriate for their description. 

We note that in some cases, if the cluster is near the border of stability, in so far as kcoui 
approaches ksurf, the vibration period increases by an order of magnitude, as for instance in 
the cases of Na4^ and Na47(j" Tabled One can expect that, under these circumstances, the 
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collective mode can be more easily excited. Experimental study of the electron relaxation and 
search for near to symmetric fission for such clusters would be of the highest interest. 

* * * 
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